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Abstract. Pulsars with drifting subpulses are thought to be an important key to unlocking the mystery of how radio pulsars 
work. We present new results from high sensitivity GMRT observations of PSR B0826 - 34 - a wide profile pulsar that exhibits 
an interesting but complicated drifting pattern. We provide a model to explain the observed subpulse drift properties of this 
pulsar, including the apparent reversals of the drift direction. In this model, PSR B0826 - 34 is close to being an aligned rotator 
Using information about the polarization and frequency evolution of the pulse profile, we solve for the emission geometry of 
this pulsar and show that the angle between the rotation axis and the dipole magnetic axis is less than 5°. As a result, our line 
of sight samples a circular path that is entirely within the emission beam. We see evidence for as many as 6 to 7 drifting bands 
in the main pulse at 318 MHz, which are all part of a circulating system of about 15 spark-associated subpulse emission beams 
that form, upon averaging, one conal ring of the mean emission. We also see evidence for a second ring of emission, which 
becomes dominant at higher frequencies (above 1 GHz) due to the nature of the emission geometry. We model the subpulse 
drift behaviour of this pulsar in detail, providing quantitative treatments of the aliasing problem and various effects of geometry 
which play an important role. The observed drift rate is an aliased version of the true drift rate which is such that a subpulse 
drifts to the location of the adjacent subpulse (or a multiple thereof) in about one pulsar period. We show that small variations, 
of the order of 3-8%, in the mean drift rate are then enough to explain the apparent reversals of drift direction seen in the data. 
We find the mean circulation time of the drift pattern to be significantly longer than the predictions of the original Ruderman 
and Sutherland (1975) model and propose an explanation for this, which relates to modified models with temperature regulated 
partial ion flow in the polar vacuum gap. The small variations in drift rate are then explained by very small heating and cooling 
effects - less than 3500 K change in the ~ 2.5 x 10* K surface temperature of the neutron star polar cap. From a detailed 
consideration of the variation of the mean subpulse separation across the main pulse window, we show that the circulating 
spark pattern is not centred around the dipole axis, but around a point much closer (within a degree or so) to the rotation axis. 
This is an indicator of the presence of a "local pole" corresponding to the non-dipolar magnetic fields that are expected to be 
present close to the neutron star surface. PSR B0826 - 34 thus provides a very rich and powerful system in which to explore 
important aspects of the physics of pulsar radio emission and neutron star magnetospheres. 

1. Introduction ical values range from 1 to 15 times Pi, the pulsar period. 

The intensity of the subpulses is systematically modulated as 

Understanding how radio pulsars work remains a challenging j^ey move along the drift band. Usually, the intensity peaks 

unsolved puzzle in astrophysics. The phenomenon of drifting j^e middle of the pulse envelope and decreases towards 

subpulses seen in some pulsars is widely thought to be an im- t^e e dges; however, in a few cases, the opposite trend is also 

portant key to understanding this mystery. This phenomenon, .gg^ (iTavlor. Manchester & Huguenin. 1975 ). The direction in 

first reported by jPrake & Craft) ^m^, manifests itself in the ^^ich the subpulses drift - either leading to trailing edge of the 

most basic form as an organised phase change of subpulses p^j^g window with pulse number, or the other way around - is 

occurring in successive pulses, resulting in the formation of another important aspect. Though both kinds of drifting pulsars 

oneor more drm bands spanning the pulse envelope jBackeij y.^.^^^ fo^ a given pulsar, the direction in which the sub- 



1973' Ranking 1986). The basic concept is illustrated in fig- puiggg jrift is usually fixed. However, there are a few excep- 

ure 8 of Backer ( 1973). In any single pulse, there are typically ^^^^^^ p^j,^., ^^ere both directions of drift are seen (though at 

two to three approximately equispaced subpulses, with a sep- different times) 
aration of P2 degrees of longitude, whose typical values are 

~ 10° - 20°. The interval between recuiTence of successive !„ the standard phenomenology (e.g. lRankinifl98^ . cases 

drift bands at a given pulse longitude is called P3 and typ- of clear subpulse driftband patterns are usually seen in pul- 

sars where the line of sight grazes along the edge of the radio 

Send offprint requests to: Y. Gupta (E-mail: ygupta@ncra.tifr.res.in) emission beam. For pulsars where the line of sight has a more 
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central traverse over the beam, the phenomenon can manifest 
itself as periodic intensity modulations of the subpulses at a 
given longitud e, without the detection of any phase changes 
llBackeilll97i . The observed periodicities related to patterns 
of drifting subpulses are usually found to be independent of ra- 
dio frequency, thus excluding all frequency dependent plasma 
effects as plausible mechanisms of the drifting subpulse phe- 
nomenon. All the properties discussed above strongly suggest 
that this phenomenon is a manifestation of the circulation of 
subpulse-associated subbeams around the magnetic axis. 

The drifting subpulse phenomenon finds a natural expla- 
nation in the model of Ruderman & Sutherland d 19751 RS75 
hereafter). This model proposes the occurrence of discrete, lo- 
calised spark discharges in a quasi-steady vacuum gap that 
forms just above the neutron star surface, as the seed activity 
that drives the radio emission mechanism. Each spark gener- 
ates a stream of electron-positron plasma that flows out along 
the dipolar field lines and, at heights that are estimated to be 
of the order of several tens of stellar radii, produces a beam of 
subpulse-associated radio emission due to some, as yet uniden- 
tified, plasma radiation mechanism. Further, this model pro- 
poses that the plasma of each spark discharge inside the vac- 
uum gap undergoes a systematic motion around the magnetic 
axis, due to the E x B drift. This, coupled with the rapidly fluc- 
tuating but rarely self-quenching nature of the sparks, leads nat- 
urally to the phenomenon of drifting subpulses. This drift rate 
is slower than the corotation rate and the sparks therefore lag 
behind with every rotation. In this picture, each drifting spark 
would correspond to one band of drifting subpulses. Since the 
spark plasma can only lag the pulsar rotation, the model pre- 
dicts drift in one direction only, depending on the relative orien- 
tation of the rotation and magnetic axes and the observer's line 
of sight for the given pulsar. The RS75 model also makes quan- 
titative predictions of the total circulation time of the drifting 
pattern for a pulsar, which, in principle, can be checked against 
observations. Further refinements of the RS75 model include 
the models of Cheng & Ruderman (1980), Gil & Sendyk (2000, 
hereafter GSOO) and Gil, Melikidze & Geppert (2003, here- 
after GMG03) which attempt to overcome some of the defi- 
ciencies and problems identified with the RS75 model. Other 
models that address the issue of subpulse drift are those o f 
Allen & Melros3 (Il983l) : Llonesl lll984l) : iKazbegi et all lll996l) : 
Wright! iI2003Ik 

However, there are some problems in the interpretation of 
the observations that make it difficult to compare with theoreti- 
cal models. Often, there are complications of the drift rates not 
being constant and stable with time. Even when the drift rates 
(or value of Pj) are stable and can be determined reliably, it is 
not clear whether they are the real values or not, as the mea- 
sured values can be modified due to the effect of aliasing. This 
effect can occur if the sampling rate (once per pulsar period) 
is too slow to track the same subpulse-associated subbeam of 
radiation in a sequence of pulses. Further, the total number of 
sparks, A^sp, needs to be known to ascertain the total circulation 
time of the pattern. In all this, the exact geometry of the ob- 
servations needs to be known for correct interpretation of the 
results. 



As a result, though there are about 40 pulsars with obser- 
vations and measurements of the phenom enon of drifting sub- 
pulses (for a summary see lRankin[ll986l) . there are only a few 
for which quantitative results have been obtained which can 
be compared against the predictions of the models. For PSR 
B0943-H10, from a careful a nalysis of the fluctuation spectra, 
lOeshpande & Rankinl ( 1 19991) were able to resolve the alias or- 
der for the observed periodicity (f 3) and also infer the presence 
of 20 spark associated subpuls e columns circulating with a time 
of 37Pi . For PSR B0809-H74. lvan Leeuwen etair(l2003h have 
recently shown that the drift rate we see is an unaliased version 
and hence P3 > 150Pi and A^sp ^ 14. For both cases, the fi- 
nal circulation time turns out to be significantly longer than the 
prediction of the RS75 model. Modified versions of this model 
(e.g. GMG03) attempt to account for these discrepancies. 

Whereas observations of pulsars with drifting subpulses 
where the line of sight just grazes a cone of emission produces 
a clear signal for analysis and interpretation, a similar phe- 
nomenon when seen for pulsars where the line of sight samples 
a large fraction of the conal ring can produce additional insight 
into the phenomenon. This makes pulsars with wide profiles 
very interesting. PSR B0826 - 34 is such a case. It has a very 
wide profile, showing emissi on over as much as 250° of lon- 
gitud e at metre wa velengths ( fPurdin et all Il979 ^. Biggs et all 
il985t iReves et all ri995.) . However, it is not a widely studied 
pulsar, proba bly because of the fa ct that it nulls for as much as 
70% of time (Durdin et al., 1979), making i t a difficult case to 
observe. The study reported bv Biggs et all ( il985.) is the most 
exhaustive study published to date for this pulsar. They found 
that there is subpulse structure at almost all longitudes in the 
wide profile of PSR B0826 - 34, with remarkable patterns of 
drifting subpulses present over 200° of longitudes. They noted 
that the drift rate shows large variations, including reversals of 
the direction of the drift - they concluded that this behaviour is 
inconsistent with the prevailing theoretical models for the phe- 
nomenon of drifting subpulses. From the integrate d pulse pro- 
files and the polarization data at 408 and 610 MHz.lBT ggs et alJ 
( 1985) concluded that PSR B0826 - 34 is hkely to be a nearly 
aligned rotator, i.e. with rotation and magnetic axes almost par- 
allel to each other, although their formal best fit values for the 
angle between the rotation and magnetic axes don't fully sup- 
port this conclusion. However, in the morphological scheme of 
iRankin (1983. 1993) it is classified as "M" type pulsar, i.e. hav- 
ing a central line of sight which goes close to the magnetic pole 
and produces distinct core and conal emission components in 
the average profile. 

It is clear from the above description that PSR B0826 - 34 
is an enigmatic pulsar that has the potential to reveal vital clues 
about the pulsar emission mechanism and therefore deserves a 
closer and more detailed study. With this in aim, we have been 
carrying out high sensitivity observations of this pulsar with 
the Giant Metre-wave Radio Telescope (GMRT). The results 
from our initial observations and the interpretations thereof are 
reported in this paper In section |2 we describe the observa- 
tions and the basic data reduction. Section|3]presents the main 
results and our model that interprets these results, including 
determination of the emission geometry for this pulsar and de- 
tailed description of various effects that modify the true values 
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of basic quantities related to the drifting phenomenon. In sec- 
tion |4] we discuss the implications of different aspects of our 
model and present a final summary. In Appendix A we sum- 
marise the basic ideas relevant to the understanding of pulsar 
emission geometry, and in Appendix B we describe the details 
of our simulation studies. 



2. Observations and Data Analysis 

PSR B0826-34 was observed on Dec 1 1, 2000 at the (GMRT), 
in the 325 MHz band, using the incoherent array mode (see 
[Oupta et al., 2000. for more details about the pulsar modes of 
operation of the GMRT). The data were obtained by incoher- 
ent addition of the dual polarization signals from 13 antennas. 
The center frequency of the observations was 318 MHz. The 
bandwidth used was 16 MHz, divided into 256 spectral chan- 
nels by the digital back-ends. The raw data were integrated to a 
time resolution of 0.5 16 milliseconds before being recorded for 
off-line analysis, where they were further integrated to obtain a 
final time resolution of 4.128 milliseconds. The duration of the 
observation was about 22 minutes, resulting in recorded data 
covering 714 periods of this pulsar. Of this, only the first 500 
or so pulses were useful; the pulsar was in the null state after 
that, and no detectable signal could be obtained from the last 
200 pulses. 

During off-line analysis, the data were first dedispersed and 
then folded synchronously with the Doppler corrected pulsar 
period to obtain the average profile, which is shown in the sec- 
ond panel of Figure IB. II The first remarkable thing about the 
average profile for PSR B0826 - 34 is the unusually large pulse 
width : emission can be seen over more than 200° of longitude. 
The "main pulse" has a double peaked structure with a compo- 
nent peak separation of 134° + 4°, centred on the longitude of 
240°. In addition, there is evidence for faint emission in the lon- 
gitude range 55°to 105°, which could be an interpulse, centred 
about 90°. Thes e features are simil ar to what is seen in the 4 08 
MHz profiles of lBiggs et all ( 1 19851) and lDurdin etlD (Il979l) . 

Comparing all the published average p rofiles available 
for this pulsar at different frequenci es (iBiggs et all 1 19851 
iDurdinetali Il979t iReves et all 1 19951 the EPN database), it 
is clear that there is a remarkable evolution of the pulse pro- 
file with frequency. It appears that the main pulse reduces in 
strength relative to the interpulse systematically, as one goes 
from lower to higher frequencies, to the extent that the inter- 
pulse actually dominates over the main pulse at frequencies 
above ~ 1 GHz (see lower 3 panels of figure IB. II for exam- 
ple). Further, the component separation between the two peaks 
of the main pulse reduces with increasing frequency, following 
the general trend seen for a large majority of pulsars (referred 
to as the "radius-to-frequency mapping" model). In Table IB. ll 
we quantify this frequency evolution by compiling all the re- 
ported values for the component separation for the main pulse 
available for this pulsar, including our new results which are 
found to be compatible with the trend. The separation between 
the centres of the main pulse and the interpulse appears to be 
roughly frequency independent, at a value of about 150° to 
160° (see figurelETl. 



As reported by iBiggs et all ( Il985l) . PSR B0826 - 34 ex- 
hibits remarkable subpulse structure in single pulse data. 
Figure lR2l shows an example of single pulse data, covering 
the first 200 pulses from our observations. Individual subpulses 
can be seen over almost the entire main pulse window. Further, 
there is evidence for regions of organised drifting of the sub- 
pulses (in the form of linear drift-bands going from later to 
earlier longitudes with increasing pulse number - we refer to 
this as negative drift) over limited period ranges. There is also 
evidence for regions of oppositely oriented drift-bands (posi- 
tive drift rate), as well as clear signs of curved drift-bands. In 
addition, there are sporadic, short duration nulls interspersed in 
the data. The over-all picture thus looks somewhat random and 
chaotic, a nd not very differe nt qualitatively from what was re- 
ported bvlBi ^gs et~all lll985ll for this pulsar at 645 MHz (their 
Fig. 2). 

However, the better signal-to-noise ratio of our data allows 
a clearer picture of the drifting subpulse pattern to emerge, 
from a gray-scale plot of the single pulse data (Figure \B3i . 
Here, one can clearly see evidence for multiple drift-bands 
across the pulse window; 6-7 bands can be seen in most of the 
sections where the pulsar does not null. However, the drift be- 
haviour of these bands of subpulses is not a regular, monotonic, 
linear drift pattern; instead it is much more irregular, exhibiting 
significant changes in the apparent drift rate, including change 
of sign. Nevertheless, the pulse regions showing drift rates of 
opposite sign are almost always connected by regions showing 
a smooth transition of the drift rate. This can clearly be seen for 
most cases where the apparent drift rate changes from negative 
to positive values, going smoothly through zero (e.g. see the 
regions of pulse # 20 to 70, pulse # 120 to 180 and pulse # 350 
to 390). For the opposite transition (from positive to negative 
drift values), the change is usually more abrupt and not always 
easily traceable. Such regions (e.g. pulse # 80 to 120 and pulse 
#180 to 220) corr e spond to the "run-away" sections identi- 
fied bvlBi ^gsetaLllll985l) . The region of pulse # 390 to 440 
is an exception, where the complete transition can be traced. 
Therefore, in the most general case, the drift pattern for this 
pulsar consists of connected regions of drift-bands of opposite 
slopes. The typical values for the maximum positive and neg- 
ative drift rates are -i-l°9/pulse (about -i-l°0/sec) to -0'!8/pulse 
(about -0°4/sec). Further, the different curved sections of drift 
bands can be connected with a continuous curve, except for 
missing sections due to nulls. A typical time scale for one cy- 
cle of drift rate reversal is about 100 periods, and about 5 such 
cycles are present in the sequence of 500 pulses in Figure Ib31 
It is also very significant that the behaviour of all the visible 
drift bands is highly correlated, i.e. all of them simultaneously 
show the same kinds of changes in drift rate. 

Since there are about 7 drift bands present in the main 
pulse region of about 160 degrees, the typical separation be- 
tween drift bands works out to about 23 degrees. A quantita- 
tive estimate of the separation between adjacent drift-bands - 
traditionally referred to as P2 - is obtained from the autocor- 
relation function of the single pulses, computed over the en- 
tire pulse window and averaged over all the available pulses. 
The result is shown in Figure IbT41 where the secondary peak is 
due to the correlation between adjacent drift-bands. From the 
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location of this peak, the average value of f 2 is found to be 
24.9 + 0.8 degrees, close to the rough estimate given above. 
This is s omewhat sma ller than the value of 29.0 + 2.0 obtained 
bv .Biggs et al.l ( fl985l) at 645 MHz. In addition, there is some 
indication from our raw data itself (e.g. Figure [B3t . of varia- 
tion of the value of f 2 across the pulse window, with a tendency 
to increase from the centre to the edge of the main pulse. We 
return to this aspect in detail in section l331 



3. Results and Interpretation 

To summarise the main results from our analysis : 
PSR B0826 - 34 has a unusually wide average profile and the 
single pulse data show multiple drifting bands of subpulses 
which drift in a highly correlated manner but show significant 
changes in the apparent drift rate, including reversals of the 
drift direction. Our basic model to explain this behaviour is 
as follows: This pulsar is an almost aligned rotator (i.e. dipole 
magnetic axis and rotation axis are almost aligned). As a result, 
our line of sight traces out a circular track through the emis- 
sion beam, centred very close to the magnetic axis. This allows 
us to sample the radiation from a large number of subpulse- 
associated plasma columns present in one particular conal ring 
of the emission beam, resulting in the large number of drift- 
ing bands of subpulses that are seen. The subpulse-associated 
plasma columns are produced by an arrangement of spark dis- 
charges occurring in the vacuum gap above the polar cap, as 
envisaged in models such as in RS75 and GSOO. The number 
of such sparks that need to be present in this particular ring, to 
explain the observed number of drift bands, is about 15. These 
spark discharges circulate in closed tracks on the polar cap, pre- 
sumably due to the E x B drift. Further, the observed drift rate 
of the these bands is not the intrinsic drift rate, but an aliased 
version. In this case, the intrinsic drift rate is such that during 
one rotation period, a subpulse drifts to reach the location of the 
adjacent subpulse, or a multiple thereof; i.e. P3 - traditionally 
defined as the time interval required for the subpulse pattern 
to repeat at the same longitudes - has a value x Pi/k (where 
k is an integer and Pi is the pulsar period). If the drift rate 
is such that f 3 = Pi/k exactly, then the observed drift-bands 
would be longitudinally stationary, i.e. zero apparent drift rate. 
Now in addition, if there are small fluctuations in the intrin- 
sic drift rate around this mean value, then the slowing down 
and speeding up of the intrinsic drift rate would be seen as ap- 
parent drifts in opposite directions. Such stroboscopic effects 
would then give rise to the apparent curved drift-bands that are 
observed in this pulsar The small fluctuations of the drift rate 
could be produced by thermal effects on t he polar cap that reg- 
ulate the ion flow in the acceleration gap jcheng & Rudermanl 
Il98(j . GMG03). 

We now discuss different aspects of our model in detail. We 
start with a determination of the emission geometry for this pul- 
sar, using information about its polarization properties and fre- 
quency evolution of the profile (section im . This is followed 
by a detailed modeling of the drift behaviour (section l3T2t . in- 
cluding the effect of aliasing and the effect of geometry on the 
estimation of f 2 - In the light of these, we present our final inter- 



pretation of the data (section]^}, deriving the possible values 
for the various drift related parameters for this pulsar. 

3.1. Determination of the geometry 

Determination of the emission geometry is vital for proper in- 
terpretation of the results. The most important parameters here 
are the following: a, the angle between the rotation and mag- 
netic axes; /?, the angle of closest approach between the mag- 
netic axis and the line of sight (also called the impact parame- 
ter); and p, the radius of the emission beam which, in general, 
is a frequency dependent quantity. The beam radius can also be 
replaced by specifying the parameters s (normalised distance 
from the dipole axis to the foot of the dipolar field line, a fre- 
quency independent quantity) and r^,,,, (the emission height of 
the radio radiation, a frequency dependent quantity). The basic 
equations describing the relationships between these quantities 
are recapitulated in Appendix A. 

Traditionally, the values of a and /3 for a pulsar are deter- 
mined from fits to the polarization angle (P.A.) curve of the 
pulsar, using the rotating vector mo del (eqnLA4|l. However, as 
has been shown by several authors jMiller & HamiltoiJ.ll993l 
IvanOmmen et all 1 19971) . there are problems in obtaining re- 
liable, independent estimates of these angles from this pro- 
cedure, in that a wide range of combination of values of a 
and /3 are found to give acceptable fits. For the specific case 
of PSR B0826 - 34, Biggs et al. ( 1985) have shown that the 
polarization angle swing is compatible with a large range of 
choices of a and Furthermore, this method does not provide 
any information about the beam radius. The latter is usually es- 
timated from a suitable measurement of the pulse width (such 
as separation between component peaks, 10% width etc) and 
the use of ean. lA.ll after the values for a and /3 have been de- 
termined (see Appendix A). 

To determine the emission geometry for PSR B0826 - 34, 
we have followed a slightly different technique which com- 
bines the available polarization information with the measured 
frequency evolution of the component separation, to constrain 
the values of the geometry parameters. Specifically, we use 
the following constraints: value of the steepest gradient of the 
P.A. curve, which we estimate to be 2.0 ± 0.5 degrees per de- 
gree ofj^ngitudejjromfigu^^ of .Biggs et al, (1985) (we note 
fliat lLvne & Manchesteil ( ll988h quote l°7/° for this value); flie 
total swing of the polarization angle over the main pulse re- 
gion, which we es timate to be a; 90°, again from figure 6 of 
iRjggsetalJ il l 9851) : the measured values of the separation of 
the main pulse components, as listed in our Table IB. ll We per- 
formed a grid search in values of a, fi and s to find all possible 
combinations of values that satisfy these constraints at all the 
frequencies listed in Table IB. II (relevant details are described 
in Appendix A). Viable solutions for a and (3 were found for 
values of the s parameter in the range 0.25 < s < 0.85. These 
solutions were in the range 1?5 < a < 5?0 and (f.6 < f3 < 2°0, 
with the value of a greater than /3 such that the requirement 
of equation iA.5\ was satisfied. The extremes of the solution 
range were : s - 0.25, a - 1°5, /3 = (f.6 for the small- 
est sized conal ring to i = 0.85, a = 5°.0, (5 - 2°.Q for the 
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largest sized conal ring. F urther, from the typical S-shaped P. A. 
curve (|Biegs et all ll985L and also the upper panel in our fig- 
ure IB.U . it is clear that we have a case of an "outer" line of 
sight, i.e. a case where the magnetic axis lies between the rota- 
tion axis and the line of sight, giving a positive value for p (see 
iNaravan & Vivekanandill982h . 

It is indeed remarkable that we are able to obtain such good 
fits matching the profile evolution over the frequency range 
of 82 MHz to 610 MHz, while meeting the basic polariza- 
tion angle requirements. Based on fits to polarization data only, 
[Biggs et al. ( 1985) obtained very different values for a and p 
for this pulsar. Using their values, it is impossible to repro- 
duce the observed pulse widths for any realistic emission height 
model. In the allowed range of 0.25 < s < 0.85, we tend to 
favour smaller values of s (w 0.25 to 0.5), as we believe that 
there is another, outer conal ring of emission for this pulsar, 
with an s value of ^ 0.7 (see section |4}- As a typical set of 
values, we use s - 0.44, a - 2°.5, /3 = l°.Q wherever specific 
values are required in the rest of this paper. For this specific set 
of values of s, a and /3, the predicted separation between the 
component peaks is 116°, 124°, 130°, 152°, 158° and 159° 
at frequencies of 610, 408, 318, 116, 95 and 82 MHz, respec- 
tively, which agree ver y well with the measureme nts listed in 
Table IB. ll We note that lLvne & Manchesteillll988h renort bes t 
fit values of a = 2?1, = 1?2 (their Table 1) and'Ra nkiiJj 19931) 
reports a = 3°,j6=l°l, both of which are compatible with our 
results. 

These results confirm our first hypothesis that we are look- 
ing at a pulsar that is very close to being an aligned rotator. 



3.2. Understanding the drift behaviour 

The basic observable quantities involved in understanding the 
subpulse drift pattern for a pulsar are illustrated in figure 8 of 
lBackej l fl97i . We start our discussion with the idealized case 
of a perfectly aligned rotator (a - 0°), in which case the line of 
sight traces out a circle on the polar cap, centred exactly around 
the magnetic axis. Further, we assume that the line of sight 
is exactly coincident with a conal ring of emission. In such a 
case, regularly spaced subpulses would be visible throughout 
the pulsar period. Later we modify our discussion to include 
realistic cases. The principal periodicities that can be measured 
are Pi, P2, Pi- The f 3 periodicity is often expressed in terms 
of a "fluctuation frequency", fj = I/P3. These can be related to 
the drift rate, D, as follows : 



/3 - D/P2 = DJ(PiP2) 



(1) 



where Pi and f 3 are in units of seconds (f^ is in Hz), P2 is in 
degrees, D is in degrees per second and Dp is in degrees per 
period. If there are A^sp equispaced sparks in a circular ring on 
the polar cap, then 



P2 = 



360° 

IvT" 



Further, the total circulation time of the pattern, i.e. the time it 
would take a drifting spark to complete one full cumulation, can 
be expressed as 



P4 



360° 
D 



(3) 



where P4 is measured in seconds. Thus, in principle, the mea- 
sured values of Pi, P2 and P3 can be used to find the drift rate 
and the total circulation time of the spark pattern, which can be 
compared with the predictions of theoretical models. 

In practice, however, the situation is complicated by the fact 
that the measured values of P2 and P3 (P™ and P™) need not 
correspond to the true values (P^ and P^); P2 can be affected 
by the drift rate and by the viewing geometry, and P3 can be 
affected due to aliasing. Also, some of the above equations are 
not valid for the general case of the non-aligned rotator. We 
first consider the effect of aUasing on P3 and then examine the 
factors affecting P2 and the effect of the viewing geometry. 

3.2.1 . Effect of aliasing on drift beliaviour 

Since the signal at any given pulse longitude is sampled only 
once per pulsar rotation period, the interpretation of any peri- 
odic signal measured at that longitude (such as the periodicity 
represented by P3) is subject to the constraints of this sampling 
rate. For sufficiently slow drift rates (such that P3 > 2Pi, or 
< 0.5/Pi), P™ is the same as P\. For P\ = 2Pi, a sub- 
pulse drifts to the longitude of the adjacent subpulse in ex- 
actly 2 periods. This corresponds to the Nyquist periodicity, 
with /n = 0.5/Pi being the Nyquist frequency. For faster drift 
rates (such that P!, < 2Pi, or > fy), an aliased version of the 
periodicity is seen. For P^ slightly less than 2Pi, P™ is slightly 
more than 2Pi and the apparent drift is of opposite sign than 
the true drift. As P^ reduces further, P™ keeps increasing cor- 
respondingly, till P3 = Pi, at which point P™ becomes infinite 
(/j = 2/n , = 0) - this corresponds to the case of longitude 
stationary subpulse patterns, with zero apparent drift. As P^ re- 
duces further (/^ > 2/n), begins to increase and a non-zero 
drift rate is again measurable; however, the direction of the drift 
is now the same as the true direction. In fact, the apparent drift 
direction reverses sign every time the true drift rate crosses a 
multiple of the Nyquist boundary (n fy), whereas the appar- 
ent drift amplitude goes through zero for every "even" crossing 
(n=even number) and through the maximum observable value 
for every "odd" crossing (n=odd number). In Appendix B we 
present results from simulations which illustrate some of these 
effects. 

Quantitatively, the above behaviour can be described as fol- 
lows : 



/' = 2kfy + i-D'f" 



(4) 



(2) 



sp 



where k - INT[(n + l)/2] and / = mod(n, 2), and n is the alias 
order such that nfy < < {n + 1)/n. This can be considered 
as a generalised form of eqn.[l] which special case obtains for 
n = 0. 

Thus, in order to estimate the true value of P3 (and hence 
determine A^sp and P4 in accordance with eqn.O the ahas order, 
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n, has to be resolved; only for n = is there a one to one 
correspondence between the measured and true quantities. In 
general, it is not always possible to infer the value of n from the 
observations. As described in section[T] it has been achieved for 
only 2 pulsars - PSR B0943+10 ( Deshpande & Rankin, 1999) 
and PSR B0809+74 ( van Leeuwen et aU 120031) - using extra 
information from specific features present in the data. 

For PSR B0826 - 34, our proposal is that we are seeing an 
aliased drift rate such that n = 2, 4, 6 ... (i.e. ^ = 1,2,3...) and 
the overall mean drift rate is such that /™ ^ 0. Thus /] » 2A:/n. 
In addition, there are small variations in /J with an amplitude 
I /™ 1^ /n, but large enough such that it crosses back and 
forth across the Nyquist boundary (« goes to n - 1, k remains 
the same, I changes sign), producing the observed pattern with 
apparent reversals of drift direction. 

The effect of aliasing also complicates the effort to resolve 
the direction of the true drift and compare it with model pre- 
dictions. In vacuum gap models, the E x B drift of the spark 
plasma in the gap is slower than the corotation speed, i.e. the 
sparks lag behind with respect to the rotation of the star (e.g. 
RS75). The direction of the measured subpulse drift that this 
translates to depends on the line of sight: for an "inner" line 
of sight (negative /?), the subpulses drift from the trailing edge 
to the leading edge of the pulse window with increasing pulse 
number, i.e. negative drift; for an "outer" line of sight (positive 
y6), the drift is from leading to trailing (positive drift). However, 
for a case where a is very small, such that the spin axis lies 
within the emission cone {a < p), there is only one direction of 
drift for any kind of line of sigh t: drift from leading to trailing 
(positive drift) - see iRuderma nl (fl976i) for similar arguments. 
Of course, the final observed direction of the drift is further 
modified by the aliasing effect as described earlier in this sec- 
tion. 

For PSR B0826 - 34, we have shown that we are dealing 
with an almost aligned rotator, with values of a that are smaller 
than the angular size of the ring of sparks fsection OTl . Further, 
we are quite certain that we are dealing with an "outer" line 
of sight geometry. Both these factors imply that the intrinsic 
drift rate should be positive (i.e. drift from leading to trailing 
in the pulse window). For an unaliased drift, the direction of 
the apparent drift will be same as that of the intrinsic drift. As 
shown earlier, the apparent drift reverses direction each time 
the amplitude of the intrinsic drift rate crosses a multiple of 
the Nyquist boundary. Now, in our model the intrinsic drift 
rate for PSR B0826 - 34 is such that we are close to an even 
Nyquist' boundary (« = 2, 4, 6...), with small variations around 
this value. When the intrinsic drift rate is slightly above the 
Nyquist boundary, the apparent drift direction should be same 
as the intrinsic. This means that the regions of positive drift 
in figure IB. 31 are regions where the drift rate is slightly faster 
than the mean rate and regions of negative drift are regions 
where the drift pattern slows down somewhat. Coupling this 



' In fact, for odd Nyquist boundaries (n=l,3,5,...) one would ob- 
serve clear "even-odd" intensity modulations like in PSR B0943+10 
(Deshpande & Rankin 1999, GSOO [their Fig. 2]) and PSR B2303+30 
(GSOO [their Fig. 3]), which are, however, not observed in PSR B0826- 
34 (see Appendix B for further discussion of aliasing effects). 



with the drift behaviour described in section |2 we conclude 
that the transition when the drift rate speeds up from less than 
the mean value to over the mean value is a relatively smoother 
and slower process than the opposite change when the drift rate 
slows down, which is more abrupt and faster. The implications 
of this are discussed further in sectional 

3.2.2. Effect of geometry on P2 

The measured value P™ does not correspond to the true value 
Pj, due to two effects. The first is the effect of the finite time 
taken by the line of sight to traverse the beam. As the line of 
sight moves from one subpulse longitude to the next one, the 
subpulse pattern drifts in longitude by a small amount; as a 
consequence, the measured value P™ is a "Doppler-shifted" 
version of the actual subpulse separation. This is illustrated 
schematically in Figure 1^3] It is straightforward to show that 
and P™ are related by 



P'2 = P™(1 - D\,/360°) , 



(5) 



where, in keeping with our convention, the sign of Z)p is posi- 
tive for drift from leading to trailing longitudes (i.e. along the 
direction of traverse of the line of sight), and negative for the 
other direction. This modified value of P2 needs to be used in 
equations ^ and |2l above. In our model for PSR B0826 - 34, 
taking the slowest drift rate (corresponding to the alias order 
of n = 2 or A: = 1) and the mean value of P™ found in section 
12 this correction amounts to ^ 7%, and increases by the same 
amount for every increment in the value of k. 

Finally, for the general case of the non-aligned rotator, the 
direct relation between P2 and the angular spacing of the sparks 
on the polar cap (eqn. ^ is no longer valid. This is because 
equal longitude intervals along the line of sight (which are mea- 
sured with respect to the rotation axis) do not map to equal an- 
gular intervals (magnetic azimuth) of the corresponding events 
on the polar cap ring (which are measured with respect to the 
magnetic axis). The general formula that relates azimuthal an- 
gles measured with respect to the magnetic axis (cr) with the 
rotational phase (0) is 



sin(cr) 



sin(a + j3)sin{(p) 
sin{Y) 



(6) 



where F, the angle between the magnetic axis and the line of 
sight at the rotational phase 0, is determined from equation lA.il 
Figure IB. 61 shows the variation of cr as a function of (p for a 
couple of sample cases of emission geometry that are relevant 
for PSR B0826 - 34. 

As a consequence of this effect, equation|2can not be used 
as it stands to compute A^sp- Instead, from a measurement of 
the longitudes of the subpulses, 0, (or of the intervals between 
them, A(/>, = P™), one first needs to correct for the effect of 
scaling due to drift (eqn.|5} and then map these corrected angles 
(or the intervals between them) to corresponding angular values 
(cr, or Acr,) on the polar cap ring. For the correct solution, the 
values cr,- will be equispaced, and this constant interval will be 



Acr = 



360° 



(7) 
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provided that the subpulse-associated subbeams circulate 
around the magnetic dipole axis. 

The conclusions from this section on understanding the 
drift behaviour can be summarised as follows. In order to esti- 
mate the total circulation time (P4), values of P3 and A^sp need 
to be estimated (eqn.|3}- The former requires knowledge of the 
alias order (eqn. 0}, and the latter requires knowledge of f 2 
(eqn. can be estimated from the variation of P™ (eqns. 

|5l |6l and IA.1> . provided the geometrical parameters (a, /3) 
are known and the drift rate is known, which in turn requires 
knowledge of f 3 and P^^. Clearly, these equations form a cou- 
pled set which can not be solved directly to obtain the values of 
A^sp, n (and hence of P3 and f 4); instead, some kind of iterative 
or grid search method is needed to find the allowed solutions. 

3.3. Interpretation of the data 

We now return to a closer inspection of our single pulse data 
to see what it can tell us about the nature of the variation of 
P2 across the main pulse window. The correlation analysis that 
was used to estimated the mean value of P™ (as discussed in 
section 13 was repeated over several narrow windows of pulse 
longitude (of about 40 degrees width), which were shifted sys- 
tematically across the entire main pulse window (in steps of 
about 4 degrees). For each case, the secondary peak of the cor- 
relation function was computed (where possible), providing an 
estimate of the value of P™ in that longitude range. Figure IB^ 
shows sample plots of the correlation functions for 5 different 
windows, exhibiting significantly different values of f ™. Figure 
IB. 81 (top panel) shows the variation of with pulse longi- 
tude, compiled from this correlation analysis. The missing data 
points correspond to cases where the correlation function did 
not show a clear secondary peak - these are regions of the pulse 
window where the drifting subpulse signal is relatively weak. 
The error bars for the estimates are not from any rigorous 
calculations or curve fitting; they are estimates by eye, based 
on the quality (level of noise and the shape of the curve around 
the secondary maxima) of the cross correlation functions. 

As seen in the top panel of figure lBTSI there appears to be a 
smooth trend for the variation of across the pulse window, 
except for a few longitudes where there is a sudden, sharp in- 
crease in the estimated value. These outliers (four are clearly 
identifiable in the figure) are not well understood and are most 
likely due to the window spanning across some pairs of very 
strong subpulses associated with some of the drift-bands; we 
ignore these four data points in the subsequent analysis. The 
figure shows that there is a measurable variation of P™ across 
the main pulse window, which can be seen more clearly in the 
trailing half of the profile. The range of variation is from a min- 
imum of about 2r.5 to a maximum of 27°, and is compatible 
with the mean value of 24° 9 that we obtained, in section|2l from 
the correlation function for the entire main pulse. The other in- 
teresting thing is that the minimum is not centred exactly in the 
middle of the on-pulse window, but appears at somewhat ear- 
lier longitudes, by about 15° or so (though the large error bars 
make it difficult to come to very firm conclusions about this). 



We now ask the question if our inferred geometry (as de- 
rived in section im can explain this observed dependence of 
P™ with longitude. The solid line in fig. IB. 61 shows the varia- 
tion of cr with for the case of the a,/3 - 2°.5, 1°0, which is 
in middle of the solution range found in section ITTI It is clear 
from this plot that the variation of cr is too fast for smaller val- 
ues of (f) and will lead to an over-correction of the smaller P™ 
values (in the middle of the pulse). This is true even after the 
measured subpulse longitudes are scaled to accommodate the 
effects of drift as described by eqn. [S] In fact, all the solutions 
found in section ITTI vield geometries where the line of sight 
is too "shallow" with respect to the emission ring, i.e. it is not 
sufficiently concentric with the emission ring. This is primarily 
because the value of a is too large, relative to that of /3. For such 
geometries, not only is the observed variation of P™ across the 
pulse window incompatible with the predictions, the observed 
intensity levels in the saddle region of the main pulse are also 
much higher than the predictions, i.e. the central bands of drift- 
ing subpulses should be much weaker than what is observed. 
This is illustrated by the dotted-line simulated profile in Fig. 1 
(see Appendix B for a description of the simulation technique), 
which is obviously inconsistent with the observed 318 MHz 
profile - the components are too narrow and the saddle is too 
weak. 

To explain the apparent disagreement between the emis- 
sion geometry and the observed detailed behaviour of the sub- 
pulse pattern, we propose the following refinement of our ba- 
sic model. Since the emission geometry constraints from pulse 
width and polarization data (section 13. U basically depend on 
the dipole field geometry at significant heights (roughly in the 
range 10-50 stellar radii), whereas the detailed variations of the 
subpulse structure are determined by events very close to the 
neutron star surface (in the "vacuum gap"), we propose that 
the centre of activity of these two are associated with differ- 
ent locations in the magnetosphere. More specifically, we pro- 
pose that the spark pattern is not centred around the dipolar 
axis of the pulsar, but instead, circulates around some other 
point on the polar cap. In models which propose higher order 
or "multi-polar" mag netic field structure very near to the neu- 
tron star surface (e.g. iGil. Melikidze fc Mitral Eooi and ref- 
erences therein), this alternate point of circulation could be the 
location of the "local magnetic pole" which, in general, could 
be offset from the dipole magnetic pole. In such a scenario, 
the generation of the spark pattern in the vacuum gap region 
is controlled by the local magnetic pole, with concentric rings 
of sparks being formed and circulating around this local pole. 
The subpulse-associated plasma columns generated by this ac- 
tivity then stream out into the magnetosphere along the dipolar 
field lines, which dominate the magnetic field structure above a 
certain height (certainly in the regions where we believe the ra- 
diation to originate). The location of the local pole (offset from 
the dipole) can be specified in terms of a new set of angles a 
and analogous to a and p (note that a' +0 - a+p needs to be 
satisfied), and a rotational angle specifying the rotational phase 
difference ^off between the fiducial plane containing the rota- 
tion axis and the dipole axis and the fiducial plane containing 
the rotation axis and the local pole axis. 
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From the variation of P™ (figure top panel) and the 
model curves in figure iB^l it is clear that in our model for 
PSR B0826 - 34, the local pole needs to be located signifi- 
cantly closer to the rotation axis that the dipole axis. Further, 
a rotational phase difference of about 15° is suggested by the 
location of the minimum of the P™ curve in the top panel of 
Fig. IB. 81 How can we make a more quantitative determination 
of the local pole geometry? We need simultaneous solutions 
for the following: values of a , /3 and (poff for the local pole, 
value of A^sp and the alias order k, such that the inferred values 
of angular separation of the sparks circulating around the local 
pole (we will use rj to refer to these angles, to differentiate from 
cr used for angles with respect to the dipole axis), will satisfy 
the corresponding version of equation^ To achieve this, a grid 
search in the values of A^^p, k, a and was carried out. For 
each possible solution point in the grid search, a few different 
choices of (poff were tried. For each combination, the f " curve 
(top panel in fig lB.8> was "corrected" to obtain a curve of Arj 
versus longitude (lower panels of fig lB.8> as follows. The mea- 
sured values at different phase locations (<pi) were converted 
to a pair of individual phase locations (0, + 0.5P™ -i- (f>off), and 
a correction for scaling due to drift was applied; these were then 
mapped to a pair of 77 values (for the specified geometry) and 
their difference converted to a A77 value for that 0,. Solutions 
that gave A;; close to 360°/A^sp across the entire pulse window 
were taken as acceptable solutions. The range of values of A^sp, 
k tried was 14, 1 to 17,4 and for a' and /3' it was to 2.2 degrees 
and 1.0 to 7.0 degrees, respectively. 

The results from this analysis are summarized in Table 
IB. 21 which lists various parameters for the different solutions. 
Acceptable solutions could be found for each choice of A^sp- 
Of these, physically meaningful solutions (see section |3 for a 
detailed explanation of this) are those for A^sp = 14, 15 and 16, 
and only these are included in Table lB.2l From column 2 of this 
table, it is seen that as the number of sparks is increased, the so- 
lution requires a higher value of alias order k (i.e. a larger value 
of drift rate). This is because the smaller f j from a larger A^sp 
requires a faster drift so that the f 2 values can be "stretched" to 
match the range of observed P" values. For every such combi- 
nation of A^sp, k values, there is small range of a ,0 values that 
result in equally acceptable solutions, in that they yield curves 
for A77 versus longitude that are very similar, and have a mean 
value, Arjg^.g, that is close to expected values A?7„p - 360°/A^sp- 
This range is a' ,0 ~ (TA^T.O to 1?2,5?2 with a' < always - 
this is in contrast to the case of the dipole solutions obtained 
in section im where a > /3. However, it is quite interesting to 
note that given the requirement a + /3 - a + /3, it is possible 
to find a combination of a and to match almost the entire 
range of dipole solutions in section lTTI Column 3 of table lB.2l 
lists a typical set of a ,/3 values. The middle panel of figure 
IB.8l shows the results for the Arj values obtained for these a , 
jS combinations for case of the three solutions, with (p„ff set to 
18°. The mean value of A77, computed from all the points across 
the pulse window (except the 4 outliers), is shown in column 5 
of table lB.2l As can be seen, these A//^,,^ values are very close 
to the expected values of 360°/A^sp (given in column 4 of the 
table). These results are very sensitive to small changes in val- 



ues of ff and , with changes as small as ~ (f.2 (in any one of 
the angles) leading to substantial deviations in the Arj curves. 
The solutions are also quite sensitive to the value of (poff, with 
small changes of the order of 10-15 degrees causing significant 
deviations from acceptable solutions. This is illustrated with 
the three sample results in the bottom panel in fig IB. 81 which 
are computed for 3 different values of (poff - The best results are 
obtained for (poff - 18°, and the values in table lB.2l are for this 
choice. 

The solutions for P3 and P4 are given in columns 6 and 7 
of table IB.2I both in units of seconds and in units of the pe- 
riod. Pi. The results for the drift rate are in column 8. The 
percentage change in the amplitude of the drift rate for the 
largest typical deviation mentioned in sectionals given in col- 
umn 9. The largest fractional change of 7.4% is needed for the 
case of A^sp,^ = 14, 1, and it reduces to 2.8% for the case of 
A^sp, k = 16,3. Thus, the changes in the mean drift rate required 
to explain the observed drift rate variations are only a few per- 
cent. 

4. Discussion and Summary 

1. On the location and arrangement of sparks in PSR 
B0826-34 : Our analysis and interpretation of the data for 
PSR B0826 - 34 reveals that we are looking at a pulsar that 
is very close to being an aligned rotator, with a value of a less 
than 5°, and more likely to be around 1° to 2°. Consequently, 
our line of sight traces out a circular track centred almost 
around the dipole magnetic axis. This not only gives the ex- 
ceptionally wide profile that we observe for this pulsar, but 
also allows us to sample the emission pattern along a ring 
around the magnetic axis. This radiation pattern reveals the 
presence of multiple drifting bands of subpulses (up to 7 are 
visible), which we interpret as produced by subpulse producing 
plasma columns associated with multiple sparks circulating in 
the ring. This is the first pulsar for which direct evidence for 
such a large number of sparks has been obtained - all previous 
cases show evidence for at most 2-3 drift bands visible within 
the pulse window. Hence our results and interpretation can be 
considered as another direct confirmation of the idea of rings 
of circulating sparks on the p olar cap, which on the average 
form c onal beams of emission dCil. Kiiak fc Seiradakis[fl99i 
iKiiak & Gil, 2002) . 

The other interesting result is that the measured variation of 
P2 - the longitude separation between the multiple drift-bands 
- across the main pulse window is incommensurate with the 
values of a and yS for the dipole geometry that we infer from 
the variation of main pulse width with radio frequency and the 
known polarization angle variation. This leads us to a new and 
interesting conclusion : that the pattern of drifting subpulses 
does not circulate around the dipole magnetic axis. The actual 
point of circulation (which we refer to as the "local pole"), is 
determined to be only ~ 3° away from the dipole axis; however, 
this small difference is enough to cause significant differences 
in the observed values of P2- 

What determines the location of this centre of circulation 
and how do we understand this in the context of the existing 
theoretical models? In our interpretation, this centre of circu- 
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lation of the drift pattern is a local pole produced by the more 
complicated magnetic field structure that could exist close to 
the neutron star surface. Though we believe that the neutron 
star magnetic field is practically dipolar at the heights at which 
the radio emission is generated in the magnetosphere, there is 
now a growing set of evidence that the magnetic field close to 
the neutron star surface (up to a few stellar radii) may have 
higher order multi-pole structure (see lUrpin &Gill2003[ for 
review). In fact, it could be a "sun-spot" like magnetic field 
(e.g. GSOO). It is worth noting that the higher magnetic field 
strength and the much smaller radii of curvature that such a 
field configuration provides are actually essential for the work- 
ing of the basic pair production mechanism i nside th e vacuum 
gap t hat drives the RS75 kind of models dOil & M elikidze. 
l2002i and references therein). In such a scenario, it is possible 
to have on the polar cap, a local pole of the magnetic field con- 
figuration whose location is not aligned with that of the pole 
of the dipolar field. In models with vacuum gaps close to the 
neutron star surface, it will be the structure of this non-dipolar 
magnetic field that will control the activities in the vacuum gap. 
In the model of GSOO, the location of this local pole is taken 
as the location of the central spark (anchored to the local pole) 
that constitutes the core beam of emission. The rest of the polar 
cap is then populated by spark discharges that are arranged in 
concentric circles around the central spark, leading to a nested 
cone structure for the emission beam. Within each circle, the 
sparks circulate due to the E x B drift in a manner similar to 
the original RS75 model. Though the higher order field domi- 
nates close to the neutron star surface, with increasing height it 
is thought to make a smooth transition to the dipolar field that 
dominates at higher altitudes. Thus we have a situation where 
the basic distribution of the seeds of emission (the sparks) on 
the polar cap is decided by the higher order surface magnetic 
field, while the actual height of the radio emission (and the evo- 
lution with frequency thereof) is dictated by the dipolar field at 
altitudes of typically about 50 stellar radii. 

Our results and inferences thus lend strong support to this 
picture of strong non-dipolar magnetic fields close to the neu- 
tron star surface. 

2. The two cone model for PSR B0826-34 : For most of 
the work presented here, we have concentrated on the emission 
properties of the main pulse region of this pulsar. However, as 
we noted in section |2] there is a weak interpulse seen in our 
318 MHz average profile, which becomes gradually stronger 
as we go to higher frequencies, and even dominates over the 
main pulse emission at frequencies higher than ~ 1 GHz (see 
fig. IB. 11 1. The natural explanation for this is emission from a 
second, outer conal ring with the geometry contriving such that 
the line of sight samples a large fraction of the inner ring at 
the lower frequencies (~ 300 MHz), while missing most of the 
outer ring. Then, as the emission beams reduce in angular size 
at the higher frequencies, the line of sight first comes closer to 
the inner ring (in the main pulse region) - at around 600 MHz 
- and then actually crosses beyond the inner ring at frequen- 
cies above 1000 MHz, leading to relatively weak emission in 
the main pulse window. This evolution naturally explains the 



increased "filling-up" of the bridge region in the main pulse 
window as we go from 318 MHz to 645 MHz. At the same 
time, the frequency evolution of the beams brings the outer 
conal ring closer to the line of sight in the interpulse region. 
This evolution is such that at around 600 MHz, significant part 
of the outer ring begins to be seen by the line of sight and at fre- 
quencies around 1 .4 GHz, the line of sight is actually grazing 
along the outer ring in the interpulse region. 

We have modeled this configuration of two conal rings 
using a simulation program that takes a given distribution of 
sparks on the polar cap and translates them along dipolar field 
lines to frequency dependent heights of emission in the magne- 
tosphere and generates the observed pulse profiles that would 
be seen along a specified line of sight (see Appendix B for 
details). The dashed lines in figure IB. II show the simulated 
profiles obtained for a model having two conal rings of emis- 
sion with frequency dependent radii determined by eqn. iA.3\ 
and centred around the local pole which is located at the point 
of the polar cap characterized by a'= 0^436, 0- 3°064 and 
'S>off - IT.3. As can be seen, the match with the observed pro- 
files is quite good, providing direct support for the model of 2 
conal rings. One important point worth stressing is that once 
the amplitude ratio of main pulse to interpulse is fixed at one 
frequency, then the spectral evolution of these profile compo- 
nents is almost entirely due to dipolar spreading of magnetic 
field lines in the emission region. 

Our two conal ring model for PSR B0826 - 34 makes clear 
predictions about how the average profile and single pulse char- 
acteristics should evolve with frequency and suitable experi- 
ments can be devised for testing these. For example, one of 
the interesting things that it predicts is that the position an- 
gle curve should be shifted with respect to the centre of the 
main pulse window, because of the 0,,// rotational phase differ- 
ence between the fiducial planes containing the local pole axis 
and the dipole axis. In figure 6 of Bi22S et al. ( 1985), there is 
some hint that the centre of the PA swing occurs at slightly later 
longitudes than the centroid of the main pulse, but the quality 
of the data is not good enough to ascertain this. Clearly, good 
quality polarimetric observations of PSR B0826 - 34 can help 
in testing this aspect of our model. Further, since our model 
predicts the interpulse emission region to arise from a different 
conal ring, we expect that the subpulse properties here could 
be different from those for the main pulse region. For example, 
depending on the number of sparks populating the outer ring, 
the measured value of P2 could be very different from the value 
we measure for the inner ring. Perhaps the most interesting fea- 
ture would be the value of drift rate and inferred values of f 3 
and P4 in the interpulse region at 600 to 1400 MHz. 

3. Pulse profile classification of PSR B0826-34 : It is clear 
from our work that there is a need for a revision of the ear- 
her classification of PSR B0826 - 34as a " M-type" pulsar (e.g. 
lRankinL[l99llLvne & Manchester, 1988), i.e. a pulsar where 
the line-of-sight and emission geometry are such that multi- 
ple cones of emission are visible along with the central core 
component. In our model, all the emission seen in the main 
pulse window originates from a single conal ring of emission 
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and that seen in the interpulse regions comes from the second, 
outer conal ring of emission; the line of sight thus does not sam- 
ple the central core region at all. Though the average profiles of 
this pulsar sometimes tend to show partial evidence for discrete 
emission components in the saddle region between the leading 
and trailing emission peaks in the main pulse region (as can be 
seen in figure IBJT i. we believe these are artefacts of the quasi- 
longitude stationary drift pattern of the multiple drift bands, 
especially when average profiles are made over relatively short 
intervals of time, from limited number of pulses. 

The classification of PSR B0826-34 as a "M-type" pro- 
file pulsar seemed strongly supported by the prominent circu- 
lar polarization, reversing sense in the central part of the pro- 
file (Biggs et al., 1985, their fig. 6). It therefore follows from 
our analysis of drifting subpulse patterns in this pulsar, that 
sense reversing circular polarization does not have to be as- 
sociated with core emission only. It can also appear in conal 
components, if the related subpulses (a s in fig.lB.3> are quasi- 
stationary in phase (for more details see iGil etallllQQl) . 

4. Inferences from the drift behaviour of PSR B0826-34 : 
We now discuss the results for the drifting behaviour of this 
pulsar in some detail. Not only do we see subpulse emission 
from several spark-associated plasma columns, we are also able 
to detect the presence of correlated drift pattern from all of 
these subpulses. The observed apparent reversals of drift di- 
rection, which are difficult to explain under any physical model 
of the subpulse drift phenomenon, are naturally explained as 
a combination of aliasing of the true drift rate coupled with 
small (x! 3 - 8%) variations of the true drift rate. Thus, in our 
interpretation, the intrinsic drift rate does not show any rever- 
sals in sign, but only small changes in amplitude, making the 
behaviour of this pulsar compatible with these models. The in- 
ferred values of both the mean drift rate and the proposed small 
variations are worthy of a detailed discussion, which we now 
present. 

As we have shown, a unique solution for the spark and 
emission geometries for PSR B0826 - 34 is not possible with 
the current data, as the alias order and A^sp can not be unam- 
biguously determined. Instead, for a given choice of the value 
of one of these, the value for the other that yields acceptable 
solutions can be found. The values for D, f 3 and P4 can then 
be found for these cases and compared with expectations from 
models. These results are summarised in Table 16^21 for three 
choices of the value of A^sp- The drift rates for most of the so- 
lutions are much lower than the values expected for the pure 
vacuum gap models (e.g. RS75, GMG03). Correspondingly, 
the values of P4 are substantially longer than the prediction. 
This result for PSR B0826 - 34 is similar to the results ob- 
tained for PSR B0943+10 (beshoande & Rankin, 'l999') and 
PSR B0809H-74 ( vanLe euwen et al., 2003). A natural expla- 
nation for this behaviour is provided by GMG03 who show 
that if the E x B drift is considered in a non-ideal vacuum gap, 
i.e. one in which there is a finite amount of charge flow (e.g. 
ICheng & Rudermanlll980l) . then the predicted value of the drift 
rate is reduced. The proposed agency for the control of the ion 
flow in the gap is thermal regulation, i.e. the ratio of the surface 



temperature to the ion critical temperature controls the rate of 
ion flow, and hence the electrical potential drop across the gap 
and hence the value for the E x B drift rate. This effect is char- 
acterised by a "shielding factor" by GMG03 who find typical 
values for it to be ~ 0.2. For our possible solutions, the values 
of the shielding factor (SF) are shown in column 10 of Table 
IB. 21 For the slowest drift rate solution, this factor is 0.36 and it 
increases to 0.94 (i.e. close to vacuum gap conditions) for the 
case of A^sp,^ - 16,3. From this it is clear that solutions for 
^ > 3 are not viable as the implied drift rate would be too fast 
- more than the maximum that is produced under pure vacuum 
gap conditions. 

In the above scenario, the small variations around the mean 
drift rate would then be due to small changes in the neutron 
star's surface temperature in the polar cap region. For exam- 
ple, to produce a 8% change in the drift rate, we find, using 
the formulation in GMG03, that only a 0.14% change in the 
surface temperature is needed. This means a change of about 
3500 K out of a mean temperature of ~ 2.5 x 10^ K. What could 
case such small and cyclic (though quasi-periodic at best) vari- 
ations in the polar cap temperature is an interesting and open 
question. One possibility is a kind of "self -regulating" mecha- 
nism of heating and cooling of the polar cap by increased and 
reduced spark activities. Whatever it maybe, from our earlier 
deductions (at the end of section l3. 2. 1> about the regions of 
slower and faster than normal drift rates, it would imply that the 
process of cooling below the mean temperature (which would 
lead to a faster than mean drift rate) is observed to be a rela- 
tively slower and smoother process than that the reverse pro- 
cess of heating (associated with slower than mean drift rates), 
which is found to be more abrupt and faster. Further, it is in- 
teresting that the observed drift rate fluctuations are found to 
be highly correlated for all the visible drift bands (section |2j, 
implying that the small temperature variations (if they are the 
cause of the drift fluctuations) occur simultaneously on a fairly 
large scale - across a large part of the polar cap. Here, detec- 
tion and comparison of subpulse drift (and its variations) in 
the interpulse region (which we believe is due a second ring 
of sparks on the polar cap) with that in the main pulse region 
should produce valuable additional insight. 

5. General implications for the average pulsar population: 
Finally, we turn to a discussion of the implications of slow vari- 
ations of drift rate for the average pulsar population. If the few 
percent variation of drift rate that we infer for PSR B0826 - 34 
is a generic phenomenon affecting a large number of pulsars, 
then it could explain some of the mysteries of non-detection 
of periodic drift variations in the fluctuation spectra analysis of 
several pulsars. Depending on the time scale of the drift rate 
fluctuations, compared to the pulsar period and the length of 
data taken for fluctuation spectra analysis, it is possible that 
drift rate variations could significantly, or even completely, 
blur the periodic signal. Thus it is possible that several pul- 
sars have drifting subpulses, but with some fluctuations in the 
drift rate, leading to reduced chances of detection of the phe- 
nomenon in a traditional fluctuation spectrum analysis. Even 
amongst the "classical" drifting pulsars (e.g. PSR B003 1 -07, 
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PSR B0943+10), variations in drift rate are a common ob- 
served feature. 

As a specific example, we would like to point out the case 
of PSR B0540+23, which almost certainly shows stroboscopic 
phenome na due to aliasing, sim ilar to PSR B0826 - 34. As re- 
ported by iNowakowskil (Il99ll) . the subpulse drift is irregular, 
without specified drift direction and with different speeds, in- 
cluding zero speed. As a result, there is no specified value of P3, 
and there is no signature in the fluctuation spectrum. Therefore, 
this pulsar is very similar to PSR 0826-34 (see Figs. 4, 5 and 
6 in Nowakowski, 1991), except that the pulse profile is much 
narrower, implying a larger inclination angle. 

6. Summary: Using results from sensitive single pulse obser- 
vations with the GMRT at 318 MHz, we have made significant 
progress in unraveling the complex drifting subpulse pattern 
of PSR B0826 - 34. In the process, we have shown that this 
pulsar is very close to being an aligned rotator. This explains 
the unusually wide profile and multiple drift bands seen for 
this pulsar. We have modeled the observed profiles and shown 
evidence for the presence of two conal rings of emission. A 
detailed treatment of effects of aliasing and geometry on the 
drifting subpulse pattern has been provided. We have given an 
explanation for the observed variability of the drift rate, includ- 
ing the apparent reversals of drift direction, within the frame- 
work of the currently available theoretical models. As a result, 
we have shown that small variations in the surface tempera- 
ture of the neutron star polar cap can explain the observed drift 
pattern. Finally, we have also demonstrated that the subpulse 
pattern shows evidence for the presence of non-dipolar mag- 
netic fields close to the neutron star surface. Clearly, there is a 
lot to be learnt from PSR B0826 - 34. 
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Appendix A: Pulsar Emission Geometry Formulae 

The g eometry of pulsar radiation is described by 
[Manche ster & Tavlor ] ll 19771 see their figure 10.4). At 
the rotational phase 0, the angle between the observer's line of 
sight and the dipole magnetic axis is given by 



cosF = cosffcosf + sinffsin^cosc 



(A.l) 



where ^ = a + (3, a is the angle between the rotation and the 
magnetic axes and /3 is the impact angle of the closest approach 
of the line of sight to the magnetic axis. 

The angle F expressed above is also the angle between the 
dipole axis and the tangent to magnetic field lines at points 
where the emission observed at the rotational phase (p origi- 



nates. Using equations for the dipolar field lines, this opening 
angle (also referred to as the beam radius) can be described as 



p = r.24 s[rUy)/RV'^P 



-1/2 



(A.2) 



where Pi is the pulsar period, rem(v) is the emission altitude 
(at which radiation at a given frequency v is generated) and 
R - 10^ cm is the neutron star radius. The mapping parameter 
0<s = d/rp<lis determined by the locus of dipolar field 
lines on the polar cap {s = at the pole and i = 1 at the polar 
cap edge), where d is the distance from the magnetic axis to the 
field hne on the polar cap corresponding to a certain detail of 
the pulse profile and r,, - 1.4 x 10"*/" 1 '^^ cm is the canonical 
polar cap radius. 

According to the generally accepted concept of radius-to- 
frequency mapping, higher frequencies are emitted at lower al- 
titudes than lower frequencies. The emission altitude can be 
described by (Kiiak & Gil, 2003, and references therein) 



(•/in _i_ Q^D ,,-0.26±0.09d0.07±0.{)3d {).37±0.05 
(40 + Vqjj^ P_i5 Pi 



(A.3) 



where R - 10^ cm, vghz is the observing frequency (in GHz), 
Pi is the pulsar period (in seconds) and P_i5 is the period 
derivative in units of 10"'^ s/s. 

The polarization position angle (PA.) curve is described by 



tariiff 



sin a sin ( 



sin^cosff - cos^cos^ 



(A.4) 



and the value of the steepest gradient of the PA. curve is 
given by 



#1 



sinff 



siny8 



(A.5) 



At a given frequency, vghz, for a chosen combination of 
a,/3 and s, the value of the beam radius, p, can be found using 
equations (IA.3I I and JA.2l i. Then, using r((f>) = p(v), the profile 
width corresponding to the emission peaks of the conal com- 
ponents, 2<p, can be found from eqn. dA.H and compared with 
observed values. For a chosen combination of a,/3, the steep- 
est gradient can be computed using eqn iA.5i and compared 
with observations; further, using the estimate of 0, the total 
swing of the polarization angle across the pulse. At//, can be 
computed from eqn ^"d compared with observed data. 

Combinations of a, (3 and s that meet all these criteria for all 
the frequencies listed in table IB. 1 1 were found by doing a grid 
search in the values of a,/3 and s. 

Appendix B: Simulations of aliasing effects 

To illustrate the effects of aliasing in the subpulse drift patterns 
of PSR B0826 - 34, we simulate the radiation of this pulsar, 
assuming that its single pulse structure reflects the circumfer- 
ential motion of 14 sparks at a distance d - Q.44rp - 33 meters 
from the pole. For the emission geometry we adopted a - 2°.5 
and /? = 1°0, and an emission altitude for the frequency of 318 
MHz was estimated from eqn. (IA.3I I. We use the simulation 
technique described in GSOO (and references therein). Each 
spark is modeled as a circular entity with a characteristic di- 
mension equal to the height h of the quasi-steady acceleration 
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gap formed above the polar cap. The adjacent sparks are sepa- 
rated from each other also by about h (spark centers are sepa- 
rated by 2h). Since the elementary pulsar radiation is relativis- 
tically beamed along dipolar magnetic field lines, the spark- 
associated intensity pattern can be transformed from the polar 
cap to the emission region and then to the observer according 
to a given geometry {a,/3, r^m). The spark-associated subbeams, 
related to subpulses observed in single pulses, perform a cir- 
cumferential motion due to the E x B drift of the spark plasma. 
The simulation program allows for circulation around either the 
magnetic dipole or local pole of surface magnetic field (assum- 
ing that this surface field is also axially symmetric on the polar 
cap), with a desired angular drift rate, D - 360° /f 4, where f 4 
is the period of spark circulation. 

For given values of P, P, a and /? the number of subpulses 
and their phases in a single pulse are determined by the value of 
Dp. We performed an experiment with time varying drift rates, 
starting with Dp - and incrementing it by 0°03 every rota- 
tional period P\ (this increment was chosen to make sure that 
cycles of gradual variations of the drift rate have duration of 
about lOOPi, as observed in real data). The idea was to find a 
sequence of about 100 pulses with varying Dp in a range ap- 
propriate to produce curved drift-bands similar to those visible 
in Fig. IB. 31 A sample result of our simulations is presented in 
Fig. IB. 91 with all the relevant information provided in the top 
panel and on both sides of the pulse window. The last column 
indicates the actual value of Dp - 0°03 n per period, where n 
is the rotational period number indicated on the vertical axis. 
Next to Dp we list values of P4/P1 = 360° /Dp. On the other 
side of the pulse window we give values of P3/P1 = f 4/(14Pi) 
and the column just before shows the fluctuation frequency, 
/3P1 ^Pi/Pi. 

What can we learn from simulated patterns presented in 
Fig. IB.9P First of all, pulses around n - I show that were there 
no drift {Dp = 0), the observer would clearly see 7 out of 14 
sparks in the form of longitude stationary subpulses. As the 
drift rate increases with the increasing number of periods, the 
subpulse drift with time varying rate becomes more and more 
apparent. However, up to about « = 100, the subpulse drift is 
relatively slow, non-aliased and proceeds from the leading to 
the trailing edge of the profile. This is the real drift direction 
and the observed drift-bands are formed by the same sparks. 
This is, however, not true in the region well above n = 100, 
where all kinds of stroboscopic aliasing effects become visi- 
ble. We have marked regions where the apparent drift-bands 
are formed by subpulses appearing at approximately the same 
phase every m-th period, where m - 5,4,3,2 and so on. It is 
worth noting that the number of apparent drift-bands is about 
7m. The drift-bands change the apparent drift direction due 
to the aliasing effect, every time crosses a multiple of the 
Nyquist frequency. 

Obviously, the region below n - 800 does not correspond 
to drifting subpulses in PSR 0826-34, because it shows alter- 
nating (even-odd), longitude stationary intensity modulations, 
which are not observed in this pulsar. It seems, however, that 
its drifting subpulse patterns are well modeled by the region be- 
tween n = 800 and n - 900, which could represent one cycle of 
multiple curved drift-bands visible in Fig. IB. 31 A clear- pattern 



of seven drift-bands is visible, moving in an aliased direction 
(from the trailing to the leading edge) in the first half of the 
cycle, and in the true direction (from the leading to the trailing 
edge) in the second half of the cycle. The drift direction change 
occurs at = I /Pi (or P3 - Pi), which is twice the con- 
ventional Nyquist frequency. At this stage the pattern advances 
exactly by one subbeam per one pulsar period Pi and the ap- 
parent drift-bands are formed by successive adjacent subbeams. 
The corresponding values of Pa/P\ - 14 and Dp - 25°7 per 
period. However, 50 pulses earlier, at the beginning of a cycle, 
Pa/P\ - 14.4 and Dp - 25° per period, while about 50 pulses 
later, at the end of a cycle, Pa/P\ - 13.3 and Dp - 2T.I per 
period. This means that the pattern speeds up along each cycle, 
increasing Dp by about 8%. 

Although we present simulations only for the case of 14 
sparks, the qualitative results are very similar for other values 
of Nsp. Results of simulations for 15 and 16 sparks, as well as 
those for 14 sparks presented in Fig. IB. 91 agree exactly with the 
drift rate solutions presented in Table lB.2l 
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Fig.B.l. Frequency evolution of the average profile of PSR B0826-34. The solid lines (bottom 3 panels) show the observed 
profiles at 3 different radio frequencies, including the result (at 318 MHz) from this paper. The dashed Unes (bottom 3 panels) 
show the best fit simulated profiles corresponding to sparks circulating around the local surface magnetic pole, while the dotted 
line corresponds to spark circulation around the dipole axis (see text for explanations). The top panel shows the position angle 
curve corresponding to dipolar geometry with a = 2°.5 and p = 1°0. Note that the inflexion point of the position angle curve lags 
the midpoint of the main pulse by about 15°. 
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Fig, B.2. Single pulse sequence of the first 200 pulses from the GMRT observations of PSR B0826-34 at 3 1 8 MHz. 
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Fig.B.3. Grey scale plot of single pulse data for 500 pulses of PSR B0826-34 from the GMRT observations at 318 MHz, with 
the average profile shown on top. The dark horizontal Unes between pulse numbers 50,60 and between 90,100 are due to radio 
frequency interference. 
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Fig. B.4. Autocorrelation results for the entire on-pulse window data in figure Ib31 The lower panel shows a detailed view around 
the secondary maxima (near the lag of 25 degrees) of the autocorrelation function. 
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Fig. B.5. Illustration of the concept of drift induced modification of P2. S i and S 2 are two sparks located at longitudes <po and (pi 
at time Tq (i.e. P2 = (pi - <po)- At Tq, the Une of sight (moving along Rws) is at <po, and sees 5 1. At a later time Ti, the line of 
sight is at (p2, where it sees spark 52- Hence = 02 - <Po- During the interval Ti - Tq, the sparks (moving along Dp) drift by an 
amount 02-01- 
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Fig. B.6. Variation of cr with (p (see eqn. |6lof text) for two different choices of geometry: a,/? - 2°.5, 1°0 (solid curve) and 
a,y6 - 0°8, 2?8 (dashed curve). The phase (p -Q° corresponds to the fiducial phase. 
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Fig.B.7. Sample results of the autocorrelation function for narrow pulse longitude windows of the data in figure Ib31 Curves for 
5 different window locations, centred at 269° (solid), 253° (long dashes), 245° (short dashes), 229° (dots) and 213° (dash-dot) 
of pulse longitude are shown to illustrate the variation of with longitude. 
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Fig. B.8. Top panel: Measured variation of with longitude. Middle panel: Inferred Arj values as a function of pulse longitude, 
for each of the solutions list in table IB. 2l the 3 horizontal solid lines are the expected values (360°/A^sp) for each case. Bottom 
panel: A77 values for 3 cases of A^sp, k - 15,2 solution, with (poff =18° (asterisks), (poff - 28° (triangles) and <poff = 8° (circles). 
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Fig.B.9. Simulations of the subpulse drift phenomenon with time varying drift rate Dp(n) - 0!03n per period, where n is the 
number of pulsar periods shown on the vertical axis (see text for explanation of other symbols). The quasi-periodic patterns of 
drifting subpulses observed in PSR B0826-34 (Fig. 3) are well modeled by the sector between n = 800 and n = 900, with Dp 
increasing from 25° to 27? 1 per period. 
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Table B.l. Frequency dependence of component separation. 



Frequency Separation References 



(MHz) (degrees) 



610 


112 


±4 


Biseset al. (1985) 


408 


128 


±4 


Bises et al. (1985): Durdin et al. (1979) 


318 


134 


±4 


This paper 


116 


151 


±4 


Turner (1980) 


95 


155 


±4 


Turner (198(y) 


82 


157 


±4 


Turner (1980) 



Table B.l. Spark geometry and drift rate solutions for 
PSR B0826 - 34 





k 


a',0 






P3 


P4 




AZ)p/Op 


SF 






n 


n 


n 


(sec) (P,) 


(sec) (Pi) 


(Vpulse) 


(%) 




14 


1 


0.8,2.8 


25.7 


25.64 


1.850 (l.OOPi) 


25.90 (M.OPi) 


25.7 


<7.4 


0.36 


15 


2 


0.6,2.2 


24.0 


24.06 


0.925 (0.50/'i) 


21.95 (7.50^1) 


48.0 


<4.0 


0.67 


16 


3 


0.8,3.2 


22.5 


22.52 


0.617 (0.33/'i) 


9.86 (5.33Pi) 


67.5 


<2.8 


0.94 



